J/^(ijj') production at the Tevatron and LHC at 0(afv 4 ) in nonrelativistic QCD 
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We present a complete evaluation for J/ip(ip') prompt production at the Tevatron and LHC at 
next-to-leading order in nonrelativistic QCD, including color-singlet, color-octet, and higher char- 
monia feeddown contributions. The short-distance coefficients of at next-to-leading order are 
found to be larger than leading order by more than an order of magnitude but with a minus sign 
at high transverse momentum p-r- Two new linear combinations of color-octet matrix elements are 
obtained from the CDF data, and used to predict J/ip production at the LHC, which agrees with 
the CMS data. The possibility of ^Sq 8 ' dominance and the J /tl> polarization puzzle are also discussed. 

PACS numbers: 12.38.Bx, 13.60.Le, 14.40.Pq 



Nearly 20 years ago, the CDF Collaboration found a 
surprisingly large production rate of ip' at high pt[J- 
To solve the large discrepancy between data and the- 
oretical predictions, the color-octet(CO) mechanism Q 
was proposed based on nonrelativistic QCD (NRQCD) 
factorization^. With the CO mechanism, QQ pairs can 
be produced at short distances in CO (^q 8 ', 3 sf\ ^Pj 8 ') 
states and subsequently evolve into physical quarkonia 
by nonpcrturbative emission of soft gluons. It can be 
verified that the partonic differential cross sections at 
leading-order (LO) in a s behave as \/p% for 3 S^\ and 

1 /p\ for ^q 8 ' and 3 Pj 8 ' , all of which decrease at high p? 
much slower than l/p% of the color-singlet (CS) state. 
The CO mechanism could give a natural explanation for 
the observed pTdistributions and large production rates 
of ip 1 and J 0|. 

However, the CO mechanism seems to encounter dif- 
ficulties in explaining the observed J/ip(ip') polariza- 
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tions. Dominated by gluon fragmentation to S[ , the LO 
NRQCD predicts transverse polarization for J/tp(ip') at 
high px |J] whereas measurements at the Fermilab Teva- 
tron give almost unpolarized J/ip(ip')\jj^. To exploit the 
underlying physics, several efforts have been made, either 
by introducing new channels @ or by proposing other 
mechanisms Q. It is a significant step to work out the 
next-to-leading order (NLO) QCD correction for the CS 
channel, which enhances the differential cross section by 
about 2 orders of magnitude at high ptQj an d changes 
the polarization from being transverse at LO into longi- 
tudinal at NLOji]. Although the CS NLO cross section 
still lies far below the experimental data, it implies that, 
compared to the a s suppression, kinematic enhancement 
at high pt is more important in the current issue. This 
observation is also supported by our recent work[Io| for 
Xc production, where we find the ratio of production rates 
°f a Xc2/ '°xei can k° dramatically altered by the NLO con- 
tribution due to change of the pr distribution from l/p% 
at LO to \jp\- at NLO in the CS P-wave channels. So we 



may conclude nothing definite until all important chan- 
nels in 1/pt expansion are presented. It means the CO 
channels "fe^ and 3 Pj 8 ^ should be considered at NLO, 
while the CS channel 3 S'| 1 ' at next-to-next-to- leading or- 
der (NNLO) in a s . Among these corrections, the com- 
plete NNLO calculation for CS is beyond the state of the 
art, and the NNLO* method is instead proposed [llj], in 
which only tree-level diagrams at this order are consid- 
ered and an infrared cutoff is imposed to control soft and 
collinear divergences, and the NNLO* contributions are 
shown to be large. However, the only 1/p^ leading contri- 
bution at NNLO in CS is given by gluon fragmentation, 
which was found [l2[ to be negligible compared to the ob- 
served J/ip(ip') production data. Other NNLO contribu- 
tions may give a \/p% term. In a complete NNLO calcu- 
lation with both real and virtual corrections, infrared and 
collinear divergences are removed and these NNLO \/p% 
contributions should be smaller than the NLO \/p% con- 
tribution due to a s suppression. Therefore, to achieve 
a good description for J/ip(ip') production a complete 
NLO calculation including both CS and CO seems to be 
necessary. 

At present, NRQCD factorization formalism with the 
CO mechanism is used to describe various processes in 
heavy quarkonium production and decay. While Jhb 
production in two-photon collisions at CERN LEP2[13| 
and photoproduction at DESY HERA[i31 are shown to 
favor the presence of CO contribution, the J/ip produc- 
tion at B factories is described well using NLO CS model 
and leaves little room for CO contributions fl5| . In order 
to further test the CO mechanism, it is necessary to study 
hadroproduction and extract CO long distance matrix el- 
ements (LDMEs) at NLO. 

In view of the importance, here we present a complete 
NLO contribution to J/ip(ip') production at the Teva- 
tron and LHC, including all important CS and CO chan- 
nels. According to the NRQCD factorization formalism, 
the inclusive cross section for direct J/ip production in 
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FIG. 1: Dependence of K factors (ratios of NLO to LO short- 
distance coefficients da) on pr in J/ip(xp') direct production 
at the Tevatron. 



FIG. 2: NLO short-distance coefficients dcrfP [ f], roda^S 1 ^], 
rider fs^], and Sum — + ridafs^] as functions 

of pr at the Tevatron, where ro = 3.9, n = -0.56 and each 
contribution is divided by der["2> ] + do"[5| ]. 



hadron-hadron collisions is expressed as 

da[pp -> J/V + X] = £ *7[(cc)„]^-^ 



(1) 



where p is either a proton or an antiproton, the indices 
i,j run over all the partonic species, and n denote the 
color, spin and angular momentum (L n ) of the interme- 
diate cc states, including 3 S\ L \ sj 8 ' and P^f 1 in the 
present issue. Compared with the iS-wave channel ob- 
tained in 0, 0, [l6| , the NLO treatment of 3 pf^ is much 
more complicated. Fortunately, using the same method 
as in , we are able to get a compact expression for the 
virtual correction, which is both time-saving and numeri- 
cally stable in the final state phase space integration. For 
technical details, we refer readers to Ref . (Io| . 

For numerical results, we choose the same parameters 
as in except that here we are restricted to y/S = 
1.96 TeV and \yj/M$i)\ < 0.6 with the Tevatron, while 
\/S = 7 TeV and < 2.4 with the LHC. 

Let us first have a glance at the overall correction be- 
haviors as presented in Fig. [T] We find the K factor 
of short-distance coefficients da for Pj channels (the 
sum over J=0,l,2 weighted with a factor of 2J+1 by spin 
symmetry in nonrclativistic limit) is large but negative at 
high pt ■ As explained in (Tp| , the negative value mainly 
originated from using the MS scheme when choosing the 
rcnormalization scheme for S-wavc spin-triplet NRQCD 
LDMEs, and does not affect the physical result. Another 

nontrivial phenomenon is that, differing from other chan- 
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nels, the K factor of S[ channel is almost independent 
of pt and not larger than 1.3. This can be understood 
since the a s correction does not bring any new kinemati- 
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cally enhanced contributions for the S[ channel, and it 
implies the expansion in a s is under control once the lead- 
ing pr (scaling as l/p T ) channel is opened up. We also 
note that K factors of all other channels are just about 1 



when pt ~ 3 GeV, which can be seen in Fig. [TJ All the 
large corrections can be attributed to the enhancement 
in 1/pt expansion. 

Since we find ^P^ channels can give a l/p T term and 
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have a large K factor, the S[ channel is no longer the 
unique source for high pt contributions. In fact, for the 
short-distance coefficients defined in Eq. ([1]) the following 
decomposition holds within an error of a few percent 



da[pp] = r a dcr^S^] + r x da[°S^} 



r3o[8h 



(2) 



where we find ro = 3.9 and r\ = —0.56 for the Tevatron, 
and ro = 4.1 and n = —0.56 for the LHC. This decom- 
position in direct J/ip(ip') production at the Tevatron is 
shown in Figj^l where each contribution is divided by 

1 [81 A 3 [81 

da[S l ] + da[S[ '] to make it easy to read. As a result, 
it is convenient to use two linearly combined LDMEs 

M m = (pj/* f s m ))+ n_ $f }) } (3) 

when comparing theoretical predictions with experimen- 
tal data for production rates at the Tevatron and LHC. 

We note that, although both (Q J ^( 3 sf ] )) and dafpf ] ] 
depend on the renormalization scheme and the factoriza- 
tion scale does not. The reason is that the de- 
pendence of (O j I^ ( 3 5*[ 8 ')) is canceled by that of r%, which 
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is originated from decomposing der[Pj } at high pt with 
all information for the dependence (here we ignore the 
contribution of 3 S'| 1 ', which decreases quickly at high pt 
in LO). So r\ should be viewed as ri(MS,/ZA) but for 
simplicity we suppress these variables in the expression. 

By fitting the pt distributions of prompt ip' and J / 
production measured at Tevatron[17j in Fig.[3]and Fig.Sl 
the CO LDMEs are determined as showing in Table U 
while the CS LDMEs arc estimated using a potential 
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FIG. 3: Transverse momentum distributions of prompt ip' 
production at the Tevatron and LHC. CDF data are taken 
from Ref. [13]. The yellow bands indicate the uncertainty due 
to CO LDMEs. 




Pt (GeV) 

FIG. 4: The same as Fig. [3] but for J lip production. The 
preliminary CMS data, taken from Ref . [18] . are compared 
with the theoretical prediction. 

model result of the wave functions at the origin^. In 
the fit we introduce a p^"* and only use experimental 
data for the region pt > p™*. In Figs. [3knd 0] and the 
following analysis, we prefer to use p™ 4 = 7 GeV. 

We find the ratio R = M^/M^ is determined to 
be as small as 0.007. Based on this fit, we may conclude 
that the direct J /ip production could be dominated by 
the <Sq channel in the chosen experimental pt region. 
To achieve this conclusion, we emphasize the following 



cut 

Pt 
GeV 


H 


(o H ) 

GeV 3 


10" 2 GeV 3 


10~ 2 GeV 3 


X 2 /d.o.f. 


7 


J/1> 
</>' 


1.16 
0.76 


0.05 ±0.02 
0.12 ±0.03 


7.4 ± 1.9 
2.0 ±0.6 


0.33 
0.56 


5 


J/rl> 

i>' 


1.16 
0.76 


0.16 ±0.05 
0.17 ±0.04 


5.2 ± 1.3 
1.1 ±0.3 


3.5 
2.2 



TABLE I: Fitted color-octet LDMEs in J/tp(ip') production 
with chosen p™*. Here ro = 3.9, ri = —0.56 are determined 
from short-distance coefficient decomposition at the Tevatron. 
Errors are due to renormalization and factorization scale de- 
pendence only. Color-singlet ( 3 S [ ^ ] ) LDMEs (O h ) are esti- 
mated using a potential model result [lij. 

points on the origination of the small R. 

(1) Wc find the fitted results are not good for data 
with pt < 7 GeV, while the data for px > 7 GeV can 
be fitted very well using the determined LDMEs for both 
J/ip and ip'. We perform a y 2 analysis for comparing 
theoretical fit with experimental data with different pjp 4 . 
Values of \ 2 /d.o.f. decrease rapidly as the cut increas- 
ing from 3 GeV to 7 GeV, and x 2 /d.o.f. becomes almost 
unchanged when p™* is larger. This may be understood 
as factorization and perturbation expansion may not be 
reliable at low pt- In Fig. 0] the curvature of observed 
cross section is positive at large pt but negative at small 
Pt, with a turning point at pt ~ 6 GeV. But the the- 
oretical curvature is positive. This implies data below 
7 GeV may not be well explained in this work (even in 
perturbativc QCD) and needs further studying. Never- 
theless, as an alternative choice, we also give the fitted 
result for p™* = 5 GeV, for which M'('^ is increased by 
a factor of 3, while the price paid is x 2 /d.o.f. increases 
from 0.33 to 3.5. The results for both p™ 4 = 7 GeV and 
p™* = 5 GeV are shown in Table H 

(2) Feed-down contributions from ip' and XcJ to J/tp 
prompt production arc properly considered. Because 
m^i and m XcJ are larger than mj/^ by only a few hun- 
dred MeV, J/ip is almost motionless in the higher char- 
monium rest frame. So pt of J/ip can be expressed 
as pt ~ p' T x (mj^/mj), where p' T and rnn are the 
transverse momentum and mass of the directly produced 
higher charmonium H . LDMEs of ip' are taken from Ta- 
ble ITJ while that of y r j are chosen with relatively smaller 
values from Ref .llOl . From experimental data in Figs. [3] 
and|U and Ref. [lOj, we see that the prompt production 
Pt distribution of J/ip is steeper than that of ip' and XcJ- 
This implies that the subtraction of more fceddown con- 
tributions will lead to a steeper J/ip direct production 
distribution and hence a smaller R. 

(3) Errors come from other sources. Varying renor- 
malization and factorization scales from tot/2 to 2mr, 
where tot = \J 4to^ ± p^ typically changes both M-j 7 ^/ 

and Mq/* by 30% (Table H}. However, the ratio R is 
almost independent of changing scales, because the de- 
pendence between two LDMEs cancels each other. Vary- 
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ing the charm quark mass m c can change the values of 
both LDMEs and R, and the dependence of R on m c is 
approximately R oc m^. Thus choosing m c = 1.5 ± 0.1 
may cause an error of 20% for R. 

So, using the Tevatron data of J /if) prompt production 
for pt > 7 GeV or even pt > 5 GeV, we find very small 
values for R, or cquivalently, M-j 7 ^ <§C (see Table 

[J). If we make a simple assumption that the smallness 
of is not due to accidental cancellation between 

(0 J /^( 3 sf ] )) and (O j ^CPq ] ))> we would have an order 
of magnitude estimate for the three LDMEs 

(0 J ^( 3 sf ] )) w {0 J ^Cpf ] ))/m 2 c < (O j/ *W)>. 
This would lead to a nontrivial result that J/ip direct 
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production is dominated by the Sq channel, hence J /if) 
is mainly unpolarized, which agrees with the polarization 
measurement [H[ . 

We have also compared our prediction for prompt J / if 
production with the CMS data in Fig.|4]and a good agree- 
ment is achieved. 

As for ip', since the difference between two LDMEs is 
not as large as that of J /if), Mf ri may be dominant at 
not too high pj-; hence, ip' may be transversely polarized 
in this region. However, it should be noted that Mf r is 

always a combination of (O^Sf 1 )) and (O^' C^ 81 )) at 
NLO; thus, whether ip' is transversely polarized at high 
p T is unclear and needs further studying. 

In summary, we calculate J/ip(ip') prompt production 
at the Tevatron and LHC at (D(aj.v 4 ), including all CS, 
CO, and feeddown contributions. A large K factor of 
P-wave CO channels at high px results in two linearly 
combined LDMEs M^ r ^' ] and , which can be 

extracted at NLO from the Tevatron data. Because of 
the steep shape of experimental J /if) prompt production 
data, we get a very small M-J 7 ^', which might indicate 
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the possibility that CO S dominates J/ip direct pro- 
duction. If this is the case, J/ip will be mainly unpolar- 
ized, which may provide a possible solution to the long- 
standing J/ip polarization puzzle. 
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Note added. Soon after this work was submitted for 
publication, a similar study appeared [20j. and for all 
color-singlet and octet channels in J/ip direct hadropro- 
duction their short-distance coefficients are consistent 
with ours. 
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